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Abstract. We study the emergence of collective dynamics in the integrable 
Hamiltonian system of two finite ensembles of coupled harmonic oscillators. After 
identification of a collective degree of freedom, the Hamiltonian is mapped onto a model 
of Caldeira-Leggett type, where the collective coordinate is coupled to an internal bath 
of phonons. In contrast to the usual Caldeira-Leggett model, the bath in the present 
case is part of the system. We derive an equation of motion for the collective coordinate 
which takes the form of a damped harmonic oscillator. We show that the distribution 
of quantum transition strengths induced by the collective mode is determined by its 
classical dynamics. 
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1. Introduction 

Many-body systems show incoherent, single-particle- motion, as well as coherent 
collective motion. Historically this phenomenon received much attention in nuclear 
physics where there is a wealth of data providing information on the coexistence of 
collective excitations, such as the Giant Dipole Resonance (GDR), and single particle 
excitations [Tj. There is also strong experimental [2] and theoretical [3] evidence that 
similar effects occur in fermionic systems different from atomic nuclei. Other examples 
for collective motion are vortex-generating rotations and oscillations in Bose-Einstein 
condensates [U [5j [6] . Furthermore collective behavior can also be present in confined 
systems such as quantum dots [TJ [8]. 

Coherent, collective motion emerges out of incoherent, single-particle motion 
whenever favored by energy conditions. Statistical analysis of spectra in nuclei indicates 
that chaotic fluctuations are due to single-particle motion, while collective motion is 
predominantly regular, for a review see Ref. [9j and more recent results in Refs p^0| [TTj . 
This generic occurrence and the coexistence of the two forms of motion pose a 
fundamental challenge. Strictly speaking, in a generic many-body system there is not an 
a priori separation of the collective motion from the single-particle dynamics. Taking 
the three-dimensional Boltzmann gas with hard-wall interactions as an example, one 
observes that the dynamics in the phase space of the system is completely chaotic [12J. 
Still, we know that the system exhibits regular collective motion in the form of sound 
waves. The deep and fascinating question in this context is therefore to understand 
from first principles how the regular motion emerges out of the full phase space chaos 

Whenever collective dynamics arises on the classical level one might expect on the 
basis of quantum-classical correspondence that this phenomena should be reflected in 
the spectral properties of the corresponding quantum many-body Hamiltonian. One way 
to probe the existence of collective excitations is to couple the system to a weak external 
periodic potential V(X) exp(iut) depending on a collective mode X. The presence of 
a collective excitation can then be usually registered as a spike at certain energies in 
the distribution of the transition strengths T(E n ) between the ground and other states 
of the system. Such a large peak can be observed, for instance, in the cross section of 
electric dipole radiation in atomic nuclei at high excitation energies, when the GDR is 
excited. On a phenomenological level one can obtain such a distribution of the transition 
strengths from a doorway-type of Hamiltonian [TJ [TJ]: 

N N 

H = ^2 E nc \n c ) (n c \ + ^ H nm \n){m\ + y"V n n e K)(n| + c.c. . (1.1) 

n c =l n,m=l n c ,n 

Here, the first term describes Nq collective states \n c ) with energies E nc , the second 
term describes the environment of single particle states \n) with H nm typically modeled 
by a random matrix. The last term models the interaction V nnc between collective and 
single-particle excitations. The collective states act as doorways into the other levels of 
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the system. A recent discussion can be found in Ref. [15]. Although successful in the 
qualitative description of collective excitation in nuclei, this model does not provide any 
explanation of the physical reasons that lead to the collective behavior. We notice that 
the collective and single-particle excitations are separated here from the start, while 
collectivity is in fact an emergent phenomena. 

Having a classical Hamiltonian whose dynamics exhibits collective motion, what can 
be stated about the distribution of the transition strengths T(E n ) for the corresponding 
quantum problem? In particular, it makes sense to ask under what conditions it is 
possible to use models like (II. II) and how the parameters there are related to the classical 
problem. It is also of considerable interest to understand the role of chaos in this 
context [THj. Unfortunately, at present we are lacking a genuine "semiclassical theory" 
for the emergence of collective excitations which would allow us to tackle the problem 
starting from the corresponding classical dynamics. The main goal of the present paper 
is to provide answers to some of the questions posed above in the framework of a 
simple integrable model of linearly coupled harmonic oscillators. The integrability of 
the system simplifies the treatment immensely. It allows for a clear identification of a 
collective coordinate X and an investigation of its dynamical evolution employing an 
analogy with the Caldeira-Leggett model [17J. After we fix the collective coordinate the 
remaining degrees of freedom are considered as a bath which is internal, not external 
as in standard models of the Caldeira-Leggett-type [HI [T9], [201 EI]- As a result, it 
turns out that the time evolution of X(t) is fully governed by the equation of motion 
for a damped harmonic oscillator of some frequency Qq determined by the parameters 
of the many-body Hamiltonian. After this we show that under certain conditions on 
the Hamiltonian of the system the averaged distribution of T(E n ) is directly connected 
to the corresponding classical problem for time evolution of X(t). In particular, the 
distribution of the transition strengths T(E n ) exhibit spikes at energies E n which are 
close to the energies E n = E + nhQ — where E is the ground state energy — of the 
collective oscillations, while the width of these spikes is controlled by the classical decay 
rate 7 of these oscillations. Even though the considered model does not involve chaotic 
features it serves as a testing ground to address the emergence of collective dynamics in a 
many-body system. Furthermore, it allows to see the effect of the absence of dynamical 
chaos on the distribution of T(E n ) and set up the ground for future investigations. 

The article is organized as follows. In Sec. [2] we introduce our model and map it 
to a Caldeira-Leggett-like system. In order to illustrate the general procedure we treat 
the special configuration of two simple coupled chains in Sec. [3J In Sec. H] we derive 
the equation of motion for the collective coordinate and obtain an expression for the 
spectral density which encodes the crucial physical properties of our model. In Sec. 
Owe investigate the distribution of transition strengths between the ground state and 
excited states and relate the result with the dynamics of collective motion. 



Collectivity in an Integrable Setting 



4 



2. Coupled chains of oscillators 

In Sec. 12.11 we define the model. After defining a collective coordinate we map the 
system onto a Caldeira-Leggett-like model in Sec. 12.21 



2.1. The model 

We consider two identical chains of one- dimensional coupled harmonic oscillators each 
consisting of N particles with positions Xj, j = 1 . . . N and momenta Pj, j — 1 . . . N as 
well as Xj and pj, respectively. They are ordered in vectors x, x, p and p. The chains 
are coupled by an interaction H int . When the coupling is "switched off" i.e., H int = 
these two chains are governed by the Hamiltonians 

#1=2^ (P> P) + ( x > W x ) > #n = ^ (P> P) + ( x ' W x ) > (2- 1 ) 

where the notation (•, •) stands for the scalar product. In the coordinate representation, 
we have 

N N 

i=i i=i 

while the potential terms describing the interactions of different particles within the 
chains can be written as 

N N 

(x, Wx) = J2 XiWijXj, (x, W x) = XiWijXj. (2.3) 

i,j=l i,j=l 

We assume that such interactions are given by a shift invariant matrix W%j = 
W / (j+n)modAf (j+n)modN which, in addition, satisfies translational symmetry condition 
Ylf=i Wij = 0. This implies that for uncoupled chains the non-interacting degrees of 
freedom are phonons. 

After introducing the coupling between the two chains the total Hamiltonian of the 
system becomes 

H = H 1 + H U + H int} (2.4) 
where the interaction term 

N N N 

Hint = ^ Kij {%i Kj) = ^ x j) 2 ^ K^jXiXj (2-5) 

i,j=l i,j=l i,j=l 

is determined by positive symmetric coupling constants Kij. In what follows, we assume 
that H(x,p,x,p) is a non-negative function. This guarantees that the motion of the 
whole system remains bounded for all times. We notice that we do not make a similar 
requirement for Hi(p,q) and Hu(p,q). 
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2.2. Mapping onto a Caldeira-Leggett-like model 

In what follows we study the dynamics of the collective coordinate X defined as the 
difference between the center of masses of two chains scaled with the factor a/ N/2 

N N 

X = ^=\^Xi =Vi r (2.6) 

To this end we map the general problem (12.41) of two coupled chains of harmonic 
oscillators to a model of Caldeira-Leggett type, where X is coupled to the "bath" 
provided by the remaining degrees of freedom. This formulation provides an intuitive 
description for the dynamics of the collective coordinate in the process of transferring 
energy from X to the bath coordinates. Such an interpretation, however, has to be 
used carefully because the energy transfer happens inside the full system and a precise 
definition of the bath depends not only on the form of the Hamiltonian (12 .4p . but also 
on the choice of the collective coordinate. 

As a first step, we introduce the new set of canonical coordinates and momenta 

N N 

Ci ^ Ai n x n , Aj n p n , (2.7) 

n=l n=l 
N N 

n=l n=l 

such that Hi and Hn become diagonal 

i=i v 7 i=i v 7 

where A^ n are the elements of the matrix A that diagonalizes W 

AWA T = — Q 2 , fi = diag(wi,...,o;j V ). (2.10) 

For the translational invariant matrix W used in this model the diagonalization matrix 
A is given by [22] 

with indices m = 2, . . . N and j — 1, . . . N. We now express the interaction part of the 
two chains in the new coordinates. For the first term in equation (12.51) we obtain 

N N N 

^ ^ ^ijip^i "Ej) ^ ^ ^ ^ kiAi n Ai m (c n C m -\- c n c m ), (2.12) 

i,j=l n,m=l i=l 



where we introduced hi = Ylj=i ^ij- Treating the second term in equation (I2.5P in an 



analogous way we obtain for the interaction part of the Hamiltonian 

H mt = ( c, K a c) + ( c, K a c) - ( c, I^c) - (c, Kf>c) , (2.13) 
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where c, c are vectors with components c n , c n and K a , K 13 are the matrices defined by 

K a = A T KA, K 13 = A T KA, (2.14) 

with Kij = Sijkj. Furthermore, after transforming the coordinates and momenta 
according to 

d n = ( -^M, d n = { -^M, (2.15) 

Vn = iXn ^ n) . Vn = iXn ^ n) (2-16) 

and defining K = K a + K@, K = K a — K 13 , the interaction term can be cast into the 
form 

H int = (d,Kd) + (d,Kd). (2.17) 
With this new set of canonical coordinates the Hamiltonian becomes 

i=i x ' i=i 

N N 

~\~ ^ K nm d n d m -\- ^ ^ K nm d n d m 

n,m=l n,m=l 

= ^( M )+(d,^ + tfd) 

+ ^ (d, ffi 2 + Kd) . (2.18) 

We notice that the collective coordinate and momentum are just 

X = d 1 , P = V i (2.19) 

and the corresponding frequency is lo\ = 0. Since X couples only to the coordinates 
d n , the part of H which depends on d, fj can be disregarded when the dynamics of the 
collective mode is considered. Consequently, the relevant part of the Hamiltonian is 
given by 

p2 N N 2 

n=2 n=2 
N 



+ £ (*«m + ^f^) d n d m . (2.20) 



m,n=2 

This already strongly resembles the Caldeira-Legget model but with the non-diagonal 
bath Hamiltonian 

N-i 2 N ~ 1 

-f^bath = £ ~~2^ ~*~ ^ B nm d n+ id m+1 , (2-21) 

n=l n,m=l 

where the elements of the (N — 1) x (N — 1) matrix B are given by 

S(n-l)(m-l) = #nm + 71,171 = 2, ... N. (2.22) 
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To cast the bath Hamiltonian into diagonal form, we introduce yet another set of the 
coordinates & = J2f=i Ujidj + i, = ^2^=1 ^jiVj+ii where U is the orthogonal matrix 
diagonalizing B, 

U T BU = — Q 2 , fi = diag(o>i,...a>jv-i). (2.23) 

With this choice of coordinates i^bath reads 

N ~i / 2 \ 

n=l ^ ' 

Now we have to perform the transformation in the part of the Hamiltonian that 
represents the interaction between the bath coordinates and the collective degree of 
freedom, 

JV-l N-1 
X Mn+i = X UnmknU = X(l, £) , (2.25) 

n=l n,m=l 

where we defined the vectors k and 1 = U T \t with the components 

iV-l 

l m = U ™kn, k n = K 1(n+1) , n = 1, . . .N - 1. (2.26) 

n=l 

Putting all the expressions together we finally arrive at the following Caldeira-Legget 
form for our model 

H' = — + KllX* + X J2 Urn + E ( a£ + J • (2 ' 27) 

m=l n=l ^ ' 

This Hamiltonian describes an effective particle moving in a harmonic potential and also 
interacting with a heat bath. We emphasize again that contrary to the Caldeira-Legget 
model the bath is part of the system and not an external configuration of particles. 
The damping of the collective motion is a result of a redistribution of energy and not 
an actual loss of energy as in models with an external bath. Furthermore, our model 
possesses only a finite number of degrees of freedom which eventually causes a return 
of energy into the collective mode. This recurrence time will be much longer, however, 
than the spreading time for sufficiently large number of particles. 



3. Chain of Oscillators with Next Neighbor Coupling as an Example 

Below we illustrate the above mapping procedure for a simple example, where the 
resulting Hamiltonian (I2.27P can be written down explicitly. We consider a system 
of two chains with next neighbor interaction coupled at one point. The Hamiltonian for 
that system reads 

N 



3=1 

,2 N 



H ^— / \ — X(j + i) mo dAr) + \Xj — £(j + i) mod 7v) J, (3.1, 

3=1 
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where Kij = (a/2)5n5ji are the coupling constants. The eigenfrequencies for a free 
chain of iV oscillators with the next-neighbor interaction as in (13.1 j) are given by 

n(k - 1) 



sin 



2N 



k 



N. 



(3.2) 



with the corresponding eigenvectors given by (12.111) . As described in section I2T21 we 
define the set of new coordinates di and consider the part of the Hamiltonian H' which 
only contains the couplings between the dj's and X = d±. Straightforward calculations 
then yield 



P a n a 
H '= 7T- + ^X 2 + -f=X(a, d + # bat]l 
2m J\ 



N 



where a = (A\2, A\n) and the bath Hamiltonian is given by 
#bath = fa, V) + ( d > (y^ 2 + « a ® a ) d } 



with f2 2 being the diagonal matrix of the eigenvalues cu^ and 
tensor product. Diagonalization of the bath leads to 

p2 N N 

2m 

with the coupling coefficients 

V2a 



(3.3) 
(3.4) 

being the ordinary 



H' 



o N AT / 2 

n=2 n=2 v 



(3.5) 



Cn(a) 



v cos 2 (^; 



iv 



E 

v fc=2 



*/ 2 N 2/7r(fe-l)' 

cos ' — — - 



f^ 2 



a; 



2\2 



E 

k=2 



V 2JV / 



(3.6) 



where the implicit equation 



4a 
iVm 



N 

E 



cos 



2 f 7r(fc-l) 
2AT 



-1 



U,'7 



(3.7) 



k=2 "1 k 

yields the eigenfrequencies luj. 



4. Dynamics of the Collective Coordinate 

We return to the general case. So far we mapped the Hamiltonian system of two 
coupled chains of harmonic oscillators to the Caldeira-Leggett model. The next step is 
to consider the time evolution of the collective mode X(t) induced by the Hamiltonian 
(I2.27p . The full quantum mechanical solution of the problem would require calculating 
the time evolution for a reduced density- matrix p rd (X) of the collective coordinate. 
While such an analysis is certainly possible along the lines of Ref. [TTJ, [23], for our 
purposes it will be sufficient to consider the most basic collective dynamical properties 
captured by the time evolution of the expectation value for the quantized collective 
observable X 

(X(t)):=Ti(pX(t)) , (4.1) 
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where p is the full density matrix. In this case the problem simplifies, since one can 
deduce the time evolution equation for (X(t)) from the corresponding equation for the 
time evolution of the quantum operator X(t) [23]. It is worthwhile to mention that, since 
H' contains only quadratic terms, the resulting equation of motion for (X(t)) coincides 
with the corresponding equation of motions for the classical observable X(t) obtained 
for the classical Hamiltonian H'. Below we give a short derivation of this equation and 
analyze its solution for certain types of initial conditions for p. 
The Heisenberg equations for our system read 

k{t) = i[H',X(t)] = ^ (4.2) 

N-l 

P(t) = 1 t[H', Pit)} = -2K U X + Unit), (4.3) 



n=l 



t(t) = kH',Ut)] = ~, (4-4) 
n m 

Kit) = ^[H',Ut)} = -rruD 2 Jn(t) + l n X(t). (4.5) 
From these equations one immediately obtains 

N-l 

mX(t) + 2K n X - Unit) = (4.6) 

n=l 

and 

mi n {t)+mu%{t)-l n X{t) = 0, n = l,...N-l. (4.7) 

We now use the representation of the momentum and coordinate operators at time zero 
in terms of creation and annihilation operators 



ln(0) = \/^- n (bn + bt), UO) = ^ibn ~ 6+). (4.8) 

With these initial conditions the solution of equation (14. 7p takes the form 

t 

+ -^- f ds sm(u n (t-s))X(s). (4.9) 
mu n J 
o 

Using this to eliminate the bath- modes from the equation (14. 8p . we obtain 

~ t oo ^ 

X(t) + —X - — [ I ' dsdu <t(Q) an(Q(t - s)) X(s) = ^ (4.10) 
mm m 



o o 

where 

N-l 



I — h — 

Fit) = l nJ ^^^ nt bn + ^ n %) (4-11) 

n=l n 



Collectivity in an Integrable Setting 



10 



is the force operator that acts on the collective coordinate and 

I 

2mw 



N-l ?2 



n=l 

is the spectral density. We further rewrite the part describing the dissipation as 

t oo t 

2 r r - , 4 /• d 




/I 
— 7 (i - s) X(s)ds, (4.13) 





where we defined the damping-kernel as 

oo 

7 (t - s ) = — [ d£u costoft - s)). (4.14) 

m J Cj 
o 

After inserting this term into equation (I4.10p we arrive at 

t 



(it 2 m J m 

o 

We now use equation (14. 15j) to obtain the evolution equation for the expectation value 
( 14.1 p of X for some class of initial states p. We assume that the initial conditions for p 
satisfies 

(1(0)) = , (P(0)> = P , (b n ) = {hi) = • (4.16) 

Here we have used the notation (A) := Tr (pA) for the expectation value of an observable 
A. Under these assumptions equation (14. 15j) yields for the expectation value of X 



dt 2 



+ nl(X(t)) + j 'ds 7 (t - s)^pt = , (4.17) 



o 



where fl^ = 2Ku/m — 7(0) and the term 7(0) is a renormalization of the potential 
resulting from the interaction between the collective mode and the bath. Equation 
(I4.17P is a classical damping equation which together with the initial conditions (14.161) 
describes the time development of the collective mode. It is straightforward to see that 
one obtains precisely the same equation for classical time evolution of X under the 
classical Hamiltonian flow induced by H' if the initial conditions are fixed as 

X(0) = 0, P(0)=P , & = 0, z/ l = 0, i = l,...JV-l. (4.18) 

We notice that the entire information on the time evolution of {X(t)) is encoded in the 
damping kernel 7. If j(t) = 7o<5(t), that is, if the system has no "memory", the above 
equation describes the damped harmonic oscillator of frequency Q with the damping 
coefficient 70. 

Since (14.171) is a linear equation, we can easily construct its solution for a general 

kernel j(t). To this end we consider a slightly different equation 

~ 00 



dt 2 



+ n 2 (X(t)}+ I ds 9(t - a) 7 (t - s ) d ( X } S » = ^ 5(t), (4.19) 
J ds m 
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with the initial conditions 

(X(-oo)) = 0, (P(-oo)} = 0. (4.20) 

at time t = — oo. Equation (14.191) describes thus the system which stays at rest for 
all times t < and then gets a "kick" at the time t = 0. After this it acquires a 
momentum Pq and continues to evolve according to equation (I4.17p . Obviously both, 
equation (14. 17j) and equation (14.191) . give the same solution for positive times. We can 
solve equation (14. 19j) employing the pair of Fourier transforms 

oo oo 

(X(t)) = J X{u) e- wt duj, X{u) = J (X(t)) e iuJt dt. (4.21) 

— oo — oo 

Applying the Fourier transformation to both sides of equation (14. 19ft we find the 
following expression 

X(uS) = TTz-o — —o t^t, (4.22) 

where 7(0;) is defined as 

j(u) := / -f(s)e tujs ds. (4.23) 







Therefore the solution of the homogeneous system becomes 

00 

(X(t)) = — / — ^ — du . (4.24) 

—00 

As one can see from equations (14.231) and (I4.24p . the dynamics of the collective mode 
is encoded in the spectral density a{uj). It is thus important to relate a{ui) to the 
interaction matrix K appearing in the original Hamiltonian (I2.20p . Recalling the 
definition (14. 12j) of a and using k = Ul we obtain 




criuj) — - 

zTtmuj \ t—' lo — uj n + le 



1 f N - 1 

— Im y 

:muj \ 

\n=l 



v L J n,n 



I-rvhuj \ lo — uj n + ie 



-ImTr 



k <g> k 



r 



(4.25) 



ul _(^ B ) 1/2 + ie 

\m J 

where 1 ® 1 T , k ® k T stands for the tensor product between 1 and 1 T (resp. k and k T ). 
The last expression can be rewritten in terms of a scalar product, 

a(u) = — Im(k, k) , (4.26) 

V 7 27rmu y u t _ (®2 + 2_K r y/ 2 + it j V ' 

where Q r , K r are (N — 1) x (N — 1) matricies obtained from Q, K by deleting the first 
row and the first column, respectively. We have now a formal expression for the spectral 
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density of our general model. Two remarks are in order. First the collective coordinate 
becomes completely decoupled from the bath if and only if k = 0. Since the components 
of k can be written as 

2 N 

fc < = 71f , (4-27) 

V j = i 

the above condition is equivalent to the requirement that the ki = Ylf=i K%i take the 
same value for all i. In particular, there is no damping if = const. We notice that 
given a splitting of the interactions: Kij — JC + 5Kij into "constant" and "fluctuating" 
parts of the interaction, only SKij contributes to k. Second, by adding the term K X 2 to 
the Hamiltonian (12. 4p one can adjust the collective frequency Cl without changing the 
spectral density a. This additional term can be incorporated into Hi, Hjj, H mt such that 
the overall structural form of H remains intact. Note that this "renormalization" results 
in a shift of the spectrum Q r of the chain Hamiltonians Hi, Hn which is compensated 
by the shift of the interaction term K r by a diagonal matrix, such that the matrix B 
(resp. a) does not change. 

The form (14.261) for the density a hinders an exact treatment for a general form 
of interaction matrix K. However, if we assume that the fluctuation part of couplings 
matrix elements are small <C m\u 2 +1 — u 2 \, we can approximate the density 

function by 

N-l , 2 

a ^ = H^ 5 { U - + 2NK/m) , (4.28) 

n=l 

where {w n } is the phononic spectrum of the noninteracting chains and the fc n 's are 
determined solely by SK^. The expression (14.281) can be interpreted to the extent 
that after introducing the interaction between the two chains the phonons acquire a 
"mass" . Assuming that k n are uniformly distributed, the behavior of o~(u) is determined 
by the spectral density of the phonon frequencies u n . In particular, in the case of 
an Ohmic law distribution for the u n this leads to c(uj) ~ ujQ(uj — — JC) at low 
frequencies. Furthermore, if K, = this in turn implies that j(t) is localized at t = 
and equation (I4.17P can be approximated by the differential equation describing time 
evolution of a harmonic oscillator with a friction. 

5. Transition strengths and collective excitation 

In the previous section, we derived an equation of motion that describes the damping 
of the collective excitation. As we mentioned already, the quantum evolution governed 
by equation (I4.17P coincides with the classical evolution of X(t) if the initial conditions 
are defined in an appropriate way. In this section we consider the problem of existence 
of quantum collective states in the spectrum of the system. One way to probe such 
collective excitations is to couple the system to an external weak periodic potential 
v(X,t) ~ A(X) cos(ujt) depending on the collective variable X. Assuming that the 
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coupling is weak, the energy absorption rate in the first order perturbation theory will 
be determined by the following spectral function 

S A (u) = £ \(0\A(X)\n)f 5 (u- ^ , (5.1) 

n=l ^ 7 

with T n = \(0\A(X)\n)\ 2 being the transition strengths between the ground state with 
energy Eq and n-th state with energy E n . The collective states can then be defined, as 
states having large transition strengths T n . Accordingly, the spectral function (15.11) keeps 
the information about the existence of collective modes in the system. Equivalently, one 
can consider the Fourier transform of Sa{^>), which is given by the time correlation of 
A(X) 

S A (t) = (0\A(X(t))A(X(0))\0). (5.2) 

On an intuitive level one might expect that the averaged transition strengths T n should 
exhibit spikes for the energies E n corresponding to collective motion. Below we show 
that under certain conditions this is indeed the case and the dynamical equation (I4.17p . 
in fact, determines the form of the time correlations <SU(i)- 

5.1. Transition strengths induced by X 

Let us first consider the case of the observable A(X) = X. We calculate the time 
correlator 

S(t) = (0\X(t)X(0)\0). (5.3) 

Since we are dealing here with a system of coupled harmonic oscillators it is useful to 
consider the set of normal coordinates {q n ,Pn) where the Hamiltonian (12 . 1 8[) becomes 
diagonal [25] . 

i=l V 7 i=l V 7 

Here cji, pi are the position and the momentum operators corresponding to (qi,Pi), with 
aj, hi being the creation and the annihilation operators, respectively. The frequencies u)i 
are the eigenfrequencies of the full system. Since the connection between old coordinates 
X, {di}, {di} and new {q{\ coordinates is given by a linear transformation, we can 
assume that 

2N 

8=1 

with some coefficients q. Substituting (15.51) into (15. 3p we obtain 
S(t) = (0\X(t)X(0)\0) 

= (0| ex.p(iHt/h)X(0)exp(-iHt/h)X(0)\0) 



£ |(0|X(0)|n)| 2 expf,^^ 



n=l 
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fc 2N ~2 

= ^E^-pHM, (5.6) 

n=l n 

where we used the relations = yjK/2mtii(a\ + hi) to calculate the transition strength 
between the ground state |0i0 2 - . - 2 tv) = |0) and excited states \nin 2 . . .u 2 n) = 
Taking then the Fourier transform of S(t) leads to 

2N ^ 2N ~ 2 

S(u>) = Y,mX(0)\n)\ 2 5(uj - u) n ) = £ - °») • ^ 

n=l n=l 

Although 5(o;) is a quantum mechanical object, we will show now that it is possible to 

relate it to the dynamics of a purely classical damped harmonic oscillator. To this end 

we consider the time evolution of the collective coordinate X under the Hamiltonian H 

with the following initial conditions: 

p 

X(Q) = — , X(0) = 0, di = 0, dM = 0, Vi > 1. (5.8) 
m 

As has been explained in the previous section, the dynamical evolution of X(t) with such 
boundary conditions is governed by equation (I4.17P for the classical damped oscillator. 
On the other hand, we can express this solution in the diagonalizing coordinates q as 
follows. The time evolution of q n (t) is given by 

q n = A n sva(u n t) . (5.9) 

where the constants A n are fixed by the initial conditions (15. 8ft : 

g„(0) = A n cu n = ^c* n . (5.10) 
m 

Accordingly, for the time evolution of X(t) we obtain 

X(t) = y2~c n q n (t) = -° Vl^LsinKt). (5.11) 

Comparing equations (15. lip and (15.71) . we see that the classical quantity X{t) and the 
imaginary part of S(t) are related via 

h 2N \~ I 2 h 
S 1 (t) := lmS(t) = -— £i|L sin ^) = ~2P X{t) - (5 - 12) 

71=1 11 

Taking the Fourier transform of S±(t) yields 

7) 2N \~ I 2 h 
= 7T- E - Un) - S(u + u n )) = ~lmX{u) (5.13) 

Am • Zu n lo 

n=l 

where X(u) is given by the righthand side of equation (I4.22p . Furthermore, comparing 
this expression with (15. 7p we recognize the connection 

2h 

S(u) = 2iO{u)S 1 (u) = —e(uj)lmX(uj), (5.14) 

"o 

where 9(u>) denotes the Heaviside step function. This can be also written explicitly as 
S(u) = ^-9(u)Im (- 2 1 , ) . (5.15) 
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It is worth noticing that this expression for S(u) can also be derived using the 
fluctuation-dissipation theorem. Suppose at a certain moment a weak time dependent 
perturbation 5H = XF ext (t) is added to the Hamiltonian (12. 4p . Under this external 
perturbation the system will be driven away from the ground state. Considering 
the linear response of the system to 5H, it follows (see e.g., [21]) that the averaged 
displacement of the collective coordinate is given by 

oo 

(X(t)) = J dt' X (t-t')FUt'), (5.16) 

— oo 

where the integration kernel is given by hx(t) = -20(t)Tm(0\X(t)X(0)\0) = -20(t)Si(t). 
On the other hand, from the previous section we know that for any force -F cxt (t) (not 
necessary weak) the evolution of (X(t)) is described by the equation: 

+ n (X(t)) +Jds 7 (t - s)^^- = ——. (5.17) 
o 

Taking the Fourier transform from both sides of this expression and comparing the result 
with the Fourier transformed equation (15.161) leads then to (15.151) . 

From equation (I5.15P we clearly see that the information on the distribution of the 
transition strengths is stored in the damping kernel 7(i) of the purely classical equation 
for the time evolution of the collective mode. One should note, however, that S(u) is 
not a smooth function but a sum of distributions with wildly fluctuating strength. It is 
easy to see, for instance, that most of the states are actually not coupled at all to the 
ground state through the operator X. Thus, in order to see a structural emergence of 
collective excitations, we need to consider a smoothened version of the spectral function 
S(u) where the average is taken over some interval [u — Auj/2,uj + Au>/2], such that 
Ao; ^> 5lj, with 5u := \uj n+ i — u n \ being the difference between two adjacent frequencies. 
We can define such a smoothened spectral function as the convolution 



S<"M := " / dS, £S ^\ 5 , (5.18) 



-oo 



where the parameter e satisfies ^> e ^> 5uj. Using then the dynamical equation (14.171) 
one obtains 

&*\u) = —0{u)Tm( w — -t^] , (5-19) 

mil \Uq — (oj — is) 1 — i{uj — ie)7 e ((jj) J 

where is the smoothened damping kernel 

oo 

%{u)= /exp((iw - e)t) j(t) dt. (5.20) 



o 



In the case when the spectral density a obeys the Ohmic law, 7 £ (co>) = 70 is constant 
and we find for the averaged S(oS) the expression 

S^{u) « —6{u) (— 2 ) . (5.21) 

mir \ (fig - uj 2 ) 2 + (cj7q) 2 / 
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Here we choose the parameter e to be small compared to 70. In the case of an 
underdamped oscillator Q > 7o/2, the above expression can be conveniently represented 
through the parameters of the corresponding classical evolution of the collective 
coordinate described by equation (14.171) . Hence we have 



X(t) = JL eX p(-7ot)sm(Q t), Q = y^-^L, 7o = 7o /2 . (5.22) 
With the parameters flo, 70 equation (15.211) takes the form 

§ M = B{u) 7~~~7T (l n\ 2 ^-2 - > ( 5 - 23 ) 

27rmfi V (w - ^o) 2 + 7o \ u + fi o) 2 + 7o / 

where we dropped the index e. In a strongly underdamped regime flo ^> 70/2 the 

transition strength distribution (I5.23j) has a maximum at the frequency uo m Q « Q of 

the collective motion, and the width of the distribution is controlled by 70, see fig. ([Q). 

On the other hand, in the overdamped regime < 7o/2 the maximum is shifted away 

from Qq and the distribution becomes very broad i.e., there are no pronounced collective 

excitations. 

5.2. Transition strengths for general couplings 

We notice that the function S(u), derived in the previous section, has only one maximum 
at a frequency near Qq. Translating this into the energy domain one concludes that the 
collective excitations show up only for the first energy level Ei = E + Q h of the 
damped harmonic oscillator, rather than for all energies E n = Eq + nQoh. This is 
directly connected with the choice of the coupling A(X) and the linear nature of our 
model, since in a harmonic oscillator the transitions induced by X only happen between 
neighboring states. Let us show that for a more general choice of the coupling A(X) 
other collective excitations show up at energies E n , n > 1 of the collective oscillator 
mode. For the sake of simplicity of exposition we will first consider the case A(X) = X 2 
and then comment on the general case. We thus consider the time correlator 

S (2) (t) := (0|X 2 (t)X 2 (0)|0) - (0|X 2 (0)|0) 2 , (5.24) 

whose Fourier transform keeps information about the transition strengths induced by 
the operator X 2 , 

00 

S( 2 \ou) := ^ / alt e^S^it) 



2vr „ 

—00 

2N 



= J2\m 2 \rn)\ 2 6L-^^) . (5.25) 

m^O ^ ' 

It is easy to show that this quantity can be expressed in terms of S(u). Indeed, 
separating the collective mode into annihilation and creation parts, 

X(t) = X + (t) +X~(t), X + (t)|0)=0, (0\X-(t)=0, (5.26) 
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and using their commutation relation leads to 

S (2) (t) = (0|X 2 (t)X 2 (0)|0) - (0|X 2 (0)|0> 2 = 2S 2 (t) . 
This immediately implies 

oo oo 

S®(w) = ^ J dt exp(iLut)S 2 (t) = 2 J S{u')S{u - J)dJ . 

— oo — oo 

Using then equation (15.141) . we obtain 

S( 2 \uj) = -8 [ S^u')^ - J)dJ . 



(5.27) 



(5.2* 



(5.29) 



If a obeys an Ohmic law and if we are in the underdamped regime, the last expression 
takes the form 

2 u 



i 



i 



V27rmn ; J V(^' + ^o) 2 + 7o - ^o) 2 + 7o 

o 

( U - .J' + fi ) 2 + 7o 2 " (w-a/-fio) 2 + 7o 2 ) ^ (5 ' 30) 
The function 5'^ 2 ^ (cij) is depicted in figure (pQ). For £7 ^ 270 (i.e., strongly underdamped 




Figure 1. The dimensionless functions (7rmOg//i) §(lj), \Trm^/ 2 /hy/2j S^(co) are 
plotted on the left-hand side (red) and the right-hand side (blue) for the parameters 
D.0 = 1, 70 = 0.1. The spikes at the bottom of the figure schematically depict the states 
which are coupled to the ground state through the operator X and X 2 , respectively. 

regime) one can clearly see a spike in the vicinity of the oscillator frequency 2Q with 
the width of the spike being twice the width of S(u) for the same parameters 70, IV 
It is straightforward to generalize the above discussion to generic observables of the 

form A(X) using the Taylor expansion 

00 

A(X) = a nX n . (5.31) 

n=0 
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After substituting this into the definition of the time correlator, and applying Wick's 
thorem to the products of X(t) we obtain 

oo 

S A (t) = (0\A(X(t))A(X)\0) = Y,PnS n (t), (5.32) 

n=0 

where (3 n are some coeficients having dimension of inverse length in power 2n. Taking 
now the Fourier transform from both sides of this expression we obtain for the spectral 
function 

1 r 00 _ 

S A (w) = y \ dt ^S A {t) = ^p n S(u)*S(u)*...*S(u\ , (5.33) 

—00 n times 

where the symbol * stands for the convolution. It is quickly seen that, in the 
underdamped regime, the n-th term of the sum (I5.33P has a maximum at the vicinity 
of nfl with a width given by wy. 



6. Conclusions 



We studied collective behavior in an integrable model consisting of two coupled chains of 
harmonic oscillators. We chose the rescaled difference of the center of mass modes of the 
chains as a collective coordinate X, and mapped our system onto a model of Caldeira- 
Leggett type. The resemblance with the well-known Caldeira-Leggett model provides 
an intuitive physical picture of the energy exchange between the collective coordinate 
and the remaining degrees of freedom playing the role of the internal bath. As a result, 
the dynamics of the collective mode is described by the damped harmonic oscillator 
equation. We then relate this dynamical equation to the problem of the existence 
of collective quantum excitations in the spectrum of the corresponding quantum 
Hamiltonian. These collective excitations are probed through the transition strengths 
induced by observables A(X), depending on the collective coordinate. As we show, for 
the dynamically underdamped regime the spikes in the distribution of the transition 
strengths appear precisely at the energies E n = E + nhfl (E = ground state energy) 
of the quantized collective harmonic oscillator, while the width of the spikes is controled 
by the damping coeficient 70 of the corresponding dynamical problem. It is worth 
mentioning that based on fluctuation-dissipation type of arguments we can extend the 
present approach to any Hamiltonian system with quadratic interactions. 

One of the important features of our model is the freedom of choice for the collective 
coordinate. Note that our definition of A in a technical sense was somewhat arbitrary. In 
principle, we could take any linear combination Y = £) i=1 [pi x i + C^i) as a collective 
coordinate, and implement the same type of mapping procedure (as in the case of A) 
onto the model of Caldeira-Legget type. We would get then precisely the same equation 
of motion for Y(t), but with a different collective frequency and damping kernel 
7(t). Not every choice for Y would be, of course, appropriate in order to regard it as a 
collective coordinate. If, for instance, the resulting dynamics becomes overdamped, no 
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clear spikes will be visible at the corresponding spectral function. On the other hand, it 
seems that there exists no "unique" choice for the collective coordinate. This means the 
parameters Qq, 7o are not intrinsic properties of the considered integrable model but 
are rather affected by the definition of the collective coordinate. It would be of a great 
interest to see whether and in what form the above "semiclassical" connection between 
the classical dynamics of a collective mode and collective excitations of the corresponding 
quantum problem can be extended to a more general class of non-integrable systems. It 
is clear that some substantial differences with an integrable case must arrise when the 
dynamics of the system becomes chaotic. 
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